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1.  INTRODUCTION 

The  United  States  Air  Force  has  recently  expressed  an  interest 
in  the  cloud  and  precipitation  microphysics  in  the  vicinity  of  the 
melting  layer  in  stratiform  clouds.  This  interest  arises  from 
observations  of  nose  cone  erosion,  and  radio  transmission  diffi¬ 
culties  in  this  area  of  the  cloud. 

The  melting  layer  itself  is  defined  by  a  "bright  band"  appearing 
on  the  RHI  (range-height  indicator)  scope  of  weather  radar.  This 
bright  band,  first  explained  by  Ryde  (1M6),  is  caused  by  a  sharp 
increase  in  the  radar  reflectivity  of  precipitation  hydrometeors  as 
they  melt.  The  acquisition  of  a  liquid  coating  causes  a  sharp  rise 
in  the  power  of  the  returned  signal.  Ci.nce  most  of  the  precipitation 
melts  at  the  same  rate,  this  causes  the  appea ranee  of  a  "bright  hand" 
on  radar.  The  band  is  narrow  in  vertical  extent,  di sappearing  when 
the  increase  in  fallspeed  of  the  melton  particles  causes  them  to 
become  less  concentrated. 

The  degree  of  radar  echo  enhancement  which  is  observed,  however, 
is  far  greater  than  that  which  can  be  explained  by  melting  alone. 
Wexler  (1952)  believes  that  the  enhancement  of  the  radar  echo  is 
caused  by  coalescence  of  ice  crystals  below  the  -3°C  level.  He 
also  indicates  that  where  little  supercooled  water  is  available, 
additional  ice  nuclei  might  come  from  shedding  of  solintors  by 
dendritic  snowflakes. 

Mason  (1?55).  Lhermitte  and  Atlas  (t?63).  and  Gunn  and  Marshall 
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(1958)  have  all  indicated  that  aggregation  of  snow  crystals  is  most, 
likely  responsible  for  the  degree  of  radar  echo  enhancement  seen 
in  the  bright  band. 

More  recently,  Lo  and  Passarelli  (Ud2)  have  developed  a  new  sam¬ 
pling  method  for  studying  the  evolution  of  snowflake  size  spectra  in 
stratiform  clouds. 

Basically,  the  method  involves  a  slow,  spiraling  descent  in  an 
instrumented  aircraft.  The  aircraft  descend:,  at  approximately  the 
same  speed  as  the  falling  snowflakes,  and  is  a  1  Lowed  to  !r'.  ft  hor¬ 
izontally  with  the  mean  wind.  Thus  the  aircraft  remains  in  a  1-D 
region  and  essentially  samples  the  same  area  of  snow  as  it  grows. 

This  type  of  flight  path  has  been  termed  the  Advent,  inr  Spiral 
Descent  (A3D). 

Snowflake  size  spectra  were  measured  continuously  during  the 
descent.  These  spectra  took  the  exponential  form 

N0  =  N0e-*D  ( i ) 

whore  N-jdD  is  the  concentration  in  the  size  range  (D,  D*-dD),  ar.d 
Nq  and  A  are  the  distribution  parameters  (intercept  and  distribution 
slope,  respectively). 

By  constructing  relationships  between  and  A  ,  Lo  and  Passa¬ 
relli  were  able  to  determine  a  pattern  of  snowflake  evolution  in¬ 
volving  depositional  growth,  aggregati onal  growth,  and  breakup  of 
aggregates. 
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It  is  their  contention  that  tho  evolution  of  the  snowflake  size 
spectrum  goes  through  three  distinct  stages.  First,  at  levels  above 
-15°C,  or  so,  growth  by  deposition  dominates.  From  approximately 
-1 to  -8°C,  or  so,  the  dominant  meehanisr  is  growth  by  aggregation. 
Below  -3*0,  or  so,  breakup  of  tho  aggregates  is  most  probably  occur¬ 
ring,  as  Nq  and  A  become  relatively  stable.  This  seers  logical  in 
light  of  the  fact  that  aggregates  larger  than,  say,  7  rom  are  only 
rarely  observed. 

We  have  found  that  a  different  anproach  to  the  analysis  of  the 
data  used  by  Lo  and  Passarelli  yields  further  support  for  the  aggre- 
gation/hreakup  mechanism  in  the  near-melting  layer  vicinity. 


2.  OBJECTIVES 


The  initial  objectives  of  this  research  were  to  continue  the 
effort  begun  by  the  author  while  serving  as  an  SCEEE/UCAF  Summer 
Faculty  Research  Fellow,  and  to  refine  the  simple  model  of  snow¬ 
flake  aggregation  developed  at  that  time  by  removing  many  of  the 
arbitrary  assumptions.  Thus,  an  accurate  picture  of  snowflake  size 
distribution  at  the  0°C  level. could  be  inferred. 

In  addition,  during  the  course  of  the  research,  it  became  evi¬ 
dent  that  it  might  also  be  possible  to  use  data  from  Mr  Force 
Geophysics  Laboratory  research  flights  to  verify  the  accuracy  of 
this  simulated  aggregation/breakup  model.  Subsequently,  a  second 
objective,  that  of  confirming  the  existence  of  an  aggregation  and, 
breakup  cycle  in  snowflake  evolution  near  the  molting  layer  was 
added. 
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3.  EXPERIMENTAL  RESULTS 

The  database  for  ti.e  research  was  provided  by  the  Air  Force 
Geophysics  Laboratory,  AFSC,  Hanscom  AF3.  FA.  The  AFGL* s  .instru- 
mentod  MC-130  cloud  physics  aircraft  made  a  number  of  flights  into 
stratiform  clouds  during  the  winter  of  1979-30.  Data  from  only 
two  of  these  flights  (Flight  30-10,  25  February;  and  Flight  30-11, 

26  February)  were  used  since  these  were  the  only  flights  which  descen¬ 
ded  to  altitudes  near  the  0°C  level.  doth  of  these  flights  were 
flown  in  the  Advecting  Spiral  Descent  pattern. 

Data  gathered  from  PFS  1-D  cloud  and  precipitation  probes  were 
used  to  determine  tho  size  distribution  o'*  snowflakes  at.  1°0  in¬ 
tervals  of  altitude  separation  as  tho  aircraft  slowly  descended 
through  the  cloud.  In  all  cases,  wo  n re  referring  to  equivalent 
melted  diameters. 


It  is  assumed  that  at  any  given  time,  tho  size  distribution  of 

snowflakes  takes  the  exponential  form 


N 


D 


N0e 


-AD 


where  N^dD  is  the  number  of  particles  in  the  size  range  (D,  D*dD), 
and  Nq  and  A  are  the  distribution  parameters  (intercept  and  distrib¬ 
ution  slope,  respectively).  To  calculate  the  total  particle  concen¬ 
tration,  N^,  at  any  given  altitude,  we  integrate  Eq.  (2)  from  D=0  to 

D=». 


■ . . ...  ..  - 
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Nt  =  J  N:Je',U’dD  (2) 

o 

The  results  of  this  integration  at  ambient  ten peraturos  f mm  -30*3 
to  - 5°C  are  shown  in  Fig.  1  for  both  fights. 

Examination  of  Fig.  1  reveals  a  r.ucber  r  f  j  ntore  at :  r-a  feat  res. 
For  flight  30-10,  the  general,  trend  of  the  total  pa rt  i  c .  o  concentra¬ 
tion  is  upward,  from  a  low  concentration  of  0.5  per  cc  at  -7i°C  to 
a  peak  of  1.7?  per  cc  at  -7°C.  For  flight  30-11,  however,  t‘o  initial 
particle  concentration  is  higher,  0.^7  per  cc,  and  actually  decrea sen, 
reaching  a  value  of  0.41  per  cc  at  the  -5°0  level. 

Both  flights  show  a  number  of  pronounced  rises  and  falls  in  the 
total  particle  concentration,  especially  at  temperatures  Warner  than 
-15°C.  Since  this  is  tho  area  in  w?".i  accreya ti on  and  breakup  are 
to  be  expected,  we  have  examined  the  d: stribition  of  cart; cl  >n  by 
size  range  as  a  percentage  of  the  total  particle  concentration. 

The  results  of  this  analysis  are  shown  in  Figs.  2a  -  2f.  n  difa 
terminate  at  the  -5  3  level  since  ha’::w  this  altitude,  accretion  of 
supercooled  droplets  is  occurring,  to  •  large  extent. 

Table  1  shows  correlation  eoeff i ■ tents  ror  the  various  s: no  ran.  «? 
of  the  two  datasets,  as  well  as  that  of  the  J,  curves.  It  .  inter¬ 
esting  to  note  that,  while  the  total  narticie  conce.utra  t :  on  curves 
are  not  well  correlated  at  all  (Correlation  Coefficient"-!. 144) ,  when 
the  data  are  compared  in  terms  of  p -rcieitagu  of  total  in  a  giver,  sice 
range,  the  correlation  coefficients  are  extremely  high,  ranging  from 
0.?49  to  0.724.  In  almost  all  cases,  increases  in  tho  percentage  of 
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particles  of  a  given  size  range  at  any  specific  temperature  occur¬ 
ring  in  Flight  >0-10,  are  seen  to  occur  at  the  same  ter  pern turo  in 
Flight  30-11.  Similar  behavior  is  seen  for  decreases  *  r.  percentage. 
There  is  no  reason  to  expect  such  a  result,  given  that  the  total 
particle  concentrations  for  the  two  flights  are  so  di  s°i  m ]  ir. 

Also  interesting  is  the  fact  that  r-oth  sots  of  data,  wV  le 
starting  out  with  very  different  size  distributions  near  the  -3')°f 
level,  seem  to  evolve  to  a  similar  size  distribution  (i.e.  their 
distribution  parameters  are  almost  equal).  This  ray  indicate  that 
the  size  distribution  of  particles  entering  the  melting  layer,  in 
terns  of  percentage  of  the  total  concentration  within  a  given  size 
range,  does  not  vary  much  from  cloud  tr  cloud,  regard! ess  of  the 
distribution  at  higher  levels  in  the  oVnd.  Perhaps  the  mechanisms 
of  aggregation  and  breakup  servo  to  establish  a  size  spoctrur  of 
particles  entering  the  molting  layer  which  is  nearly  constant 
(percentage-wise)  from  one  cloud  to  the  next. 

A  comparison  of  Figs.  2a  and  2c  illustrates  the  most  striking 
example  of  the  effect  of  aggregation  .and  breakup  on  size  distribution. 
It  is  observed  that  increases  in  the  percentage  of  particles  cf 
diameter  2-3  nan  are  always  accompanied  by  decreases  in  the  percentage 
of  particles  smaller  than  1  mm  in  diameter.  Indeed,  the  curves  for 
these  two  size  ranges  show  an  almost  perfect  negative  correlation 
for  both  flights. 

It  would  appear  that  aggregation  o'  snowflakes  bogies  with  par¬ 
ticles  lass  than  1  mm  in  diameter,  proceeds  rapidly  through  the  1-2 
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Fin  diamoter  range  and  reaches  the  lower  end  nf  the  2-3  r>*  diameter 
range  with  no  significant  breakup.  lb  w ever,  breakup  appears  to  :n- 
crease  sharply  as  the  particles  appro**  :h  3  n«  in  bare*.  •••r.  n’-onrur 
breakup  occurs,  there  is  a  sharp  drop  in  the  j.-orcentar  ••  i  f  1-3  r« 
diamoter  particles  and  an  accompany: nr  rise  ir.  the  percentage  of 
particles  loss  than  1  h.f  in  diameter.  Lt  is  interesting  to  note  teat 
in  Flight  30-10,  breakup  of  2-3  r.r,  diameter  'an!  larger)  particles 
did  not  result  in  significant  increases  of  1-2  diameter  particles; 
however,  this  was  not  the  case  in  Flight  30-11,  where  increases  in 
1-2.  mm  diameter  particles  were  seen  during  breakup  of  larger  parti¬ 
cles.  This  may  be  an  effect  rebated  to  t.ho  *<  piati  z*,  tion  nrocess 
for  the  distribution  as  we  near  the  melting  layer.  indeed,  the  per¬ 
centage  of  1-2  mm  diamotor  particles  in  F’ir1*'  iO-10  r  >:•  a  ;  r.ed  essential¬ 

ly  constant  at  ZJ-Z'i'/i  from  the  -1  l-n.-el  down  t.o  -5°0,  while  the 
percentage  of  1-2  mm.  diameter  particles  in  Fli.-ht  30-1 1  did  not  be¬ 
come  stable  at  this  23 -2*£  level  urst*  1  -1  )°1. 

Larger  size  ranges  exhibited  sis':  >r  >••»*.*»  vlor  (Fir,.  2d-.*’.-,  w'  th 
periods  of  aggregation  indicated  by  upswings  in  the  percentage  of  ' arye 
particles  accompanied  by  sharp  decross  is  in  tho  nur  her  of  era  11  par¬ 
ticles,  and  periods  of  breakup  shewing  just  the  opposite  behavior. 

One  of  the  initial  problems  encountered  at  the  start  of  the  re¬ 
search  was  that  we  had  no  idea  what  size  the  fragments  of  a  broken-up 
snowflake  would  be.  Our  analysis  of  tho  data  from  those  two  flights 
seems  to  indicate  that  breakup  of  snowflakes  will  become  a  -'oworful 
mechanism  when  the  flakes  approach  3  mm  in  diameter.  The  breakups 
produce  particles  mainly  smaller  than  l  mm  in  diameter,  with  sore 
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1-2  mm  diameter  particles  produced  as  well.  Foot  snowflakes 
reaching  the  critical  size  will  brook  up  (only  ?-ll't  of  the  par¬ 
ticles  crossing  the  -5°C  level  ware  4  mm  in  d'  meter  or  larger. 
In  Flight  30-11,  only  "}$  of  the  partial  jo  entering  the  molting 
layer  were  4  mm  in  diameter  or  larger). 
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4.  COMPARISON  OF  ASD  DATA  WITH  THEORETICAL  MODEL 

A  simple  model  of  snowflake  aggregation  and  breakup  was  de¬ 
veloped  prior  to  this  research  (Newman,  1931).  This  model  con¬ 
tained  a  number  of  assumptions  which  have  now  neon  modified. 

The  basis  for  the  model  is  the  equation  for  snowflake  aggre¬ 
gation  given  by  Rogers  (1979): 

dm/dt  =  SMTtR2A-u  (4) 

where  dm/dt  is  the  mass  growth  rate,  2  the  mean  collection  efficiency, 
M  the  cloud  ice  content,  R  the  radius  of  the  collecting  snowflake, 
and  au  the  difference  in  fallspeed  between  the  collootinr  snow¬ 
flakes  and  the  collected  ice  crystals.  Mr  son  (19?l)  has  related  the 
diameter,  D,  of  stellar  dendritic  snowflakes  to  their  mass  by, 

m  =  0.027D2  (5) 

and,  substituting  into  Eq.  (4)  yields, 

l/d  dr/dt  =  4.63HMkau  {(:,) 

The  initial  model  assumed  a  constant  value  nor  Au  based  on  earlier 
research  by  Magono  (1953)  which  indicated  that  snowflake  fallspeed 
was  essentially  independent  of  sizo  for  dendritic  flakes.  However, 
Langleben  (195+)  Has  shown  that  snowflake  fallspeed  is  a  function  of 
size  such  that, 
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U  =  193.353H0,31  (?) 

We  have  further  assumed  that  the  collect')!  ice  crystals  fall  with 
constant  speed  of  30  cm/ sec. 

f)#  a] 

Thus  we  obtain  the  expression  (193.353H  *"  -39)  to  be  substi¬ 

tuted  for  au,  yielding, 

dR/R  =  4.63EMH  (191.353R0*31  -30)dt  (3) 

If  we  assume  that  M  is  relatively  constant,  and  that  the  expression 
in  parentheses  is  constant  over  a  sra] 1  interval  of  time,  t,  then 
the  final  aggregation  equation  will  be, 

R  =  Rqsxp  [4.63»;n(W3.353i^'yl-30)tj  (9) 

This  equation  can  be  solved  for  intervals  of  t=19  seconds.  All 
that  is  required  are  values  of  E  and  M. 

In  the  1931  model  runs,  E  values  of  0.1,  0,9  and  1.9  were  used. 
However,  Passarelli  (1973)  has  shown  that  there  is  a  strong  wake 
effect  associated  with  falling  dendritic  snowflakes.  This  wake  ef¬ 
fect  makes  possible,  values  of  S  greater  than  unity,  i’assaroili 
found  a  value  of  1.410.6  to  be  in  good  agreement  with  observed 
rates.  We  have  used  this  value  of  1,4  >n  the  model 

runs  here. 

In  addition  to  the  modifications  in  An,  we  have  used  the 
moist  adiabatic  lapse  rate  of  -5.33°C/km  in  these  runs.  This  is 
very  close  to  the  observed  lapse  rate  for  both  AP’GL  flights. 

The  earlier  model  had  merely  assumod  lapse  ratos  of  -4,  -5,  and 
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-6#C/km.  Thus  rany  of  the  simplifying  assumpti ons  made  in  the  ear¬ 
lier  work  have  been  modified  or  elii'instod  in  favor  of  rore  recent 

observed  data. 

Table  2  shows  the  results  of  the  model  runs.  Values  of  M  were 

3  3 

allowed  to  vary  from  0.5  g /nr  to  2.0  g/m  .  Note  that  for  a  cloud 
ice  content  of  1.5  g/m  ,  the  model  predicts  3  occurrences  of  break¬ 
up  at  3  nan  diameter  between  -15*C  and  -5°C.  These  predicted  break¬ 
ups  occur  at  -11.8°C,  -3.6  C,  and  -5«4*C.  These  temperatures  com¬ 
pare  quite  favorably  with  both  AFGL  flights,  where  breakup  is  appar¬ 
ently  occurring  at  -13°C,  -10°C,  and  -5°G.  If  the  critical  diameter 
is  assumed  to  be  as  high  as  4  mm,  3  breakups  will  occur  at  -12.0*C, 
and  -6.1°C  for  a  cloud  ice  content  of  2.0  g/nA  Again, 
these  predicted  breakup  temperatures  compare  favorably  with  observed 
data. 
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5.  SUGARY  AND  CONCLUSIONS 

Lo  and  Passarelli  (1932)  have  conducted  experiment 5  in  strati¬ 
form  clouds,  dealing  with  snowflake  size  distributions.  They  have 
used  a  sampling  method  (the  Advecting  Spiral  Descent)  which  follows 
the  snowflakes  down  through  the  cloud  at  their  aporoximate  fallspeed. 

They  have  defined  three  stages  in  the  evolution  of  snow  size 
spectra: 

1)  Depositional  growth— occurring  predominantly  at  temper- 
atures  colder  than  -13  to  -lo  C. 

2)  Growth  by  aggregation — occurring  nredominamtly  at 
temporatures  between  -1^  and  -3*C. 

3)  Breakup  of  aggregates — which  begins  at  about  -3°C 
and  continues  down  to  the  melting  layer. 

We  have  approached  this  problem  in  a  somewhat  different  man¬ 
ner,  3y  examining  the  total  particle  nun  nor  concentration,  and  the 
percentages  of  this  total  of  various  size  particles  from  less  than 
1  mm  to  more  than  5  *nm  in  diameter,  we  have  been  able  to  determine 
the  altitudes  at  which  aggregation  and  breakup  are  occurring. 

Our  analysis  of  the  data  clearly  shows  that  when  aggregation  is 
occurring,  the  percentage  of  the  total  distribution  of  large  parti¬ 
cles  increases  at  the  expense  of  the  smaller  particles.  Conversely, 
when  breakup  occurs,  there  is  a  sharp  decrease  in  large  particles, 
which  is  accompanied  by  a  rise  in  the  percentage  of  smaller  particles. 
These  fluctuations  in  size  distribution  are  well  explained  by  an 


aggregation/breakup  mechanism  such  as  that  proposed  by  Lo  and 

Passarelli. 

In  addition,  comparison  of  the  a  'parent  temperatures  at  which 
breakup  is  occurring  in  the  cloud,  to  predicted  temperatures  from 
a  simple  model  of  aggregation  and  breakup  shows  a  high  degree  of 
correlation.  This  would  appear  to  provide  further  evidence  for  the 
existence  of  such  a  mechanism. 

The  remaining  question  is — does  this  mechanism  actually  occur? 
Unfortunately,  there  are  not  enough  flight  data  available  to  provide 
a  concrete  answer.  Further  research  will  be  needed  to  verify  the 
existence  of  this  mechanism  in  clouds,  near  the  melting  layer. 
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TABLE  1 


Size  Range 

All  sizes 

Porcent  1  mm  diameter 
Percent  1-2  mm  diameter 
Percent  2-3  mm  diameter 
Percent  3-^  mm  diameter 
Percent  4-5  mm  diameter 
Percent  5  mm  diameter 


Correlation  Coefficient 

-0.1440 

40.9439 

+0.3336 

+0.935? 

+0.3318 

+0.3244 

*•0.7233 


Table  1.  Correlation  coefficients  for  various  size  ranfte^  Vietwcen 
datasets  30-10  -and  80-11 
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Critical 
_ £mi 


3 


4 


Table  2. 


TABLE  2 


diameter  Cloud  ice  content  Breakup  terperature 

l) _  (g/™3)  _ (  C)  _ 


0.5 

-5.5 

1.0 

-10.3,  -5.5 

1.5 

-11. 3, -3. 6, -5.4 

2.0 

-12. 6, -10. 2,-7. 3t- 

0.5 

1.0 

-9.1 

1.5 

-11.1, -7.1 

2.0 

-12.0, -9.0, -6.1 

Predicted  breakup  temperatures  Tor  critical  diameters  of 
3  and  4  mm  at  various  cloud  ice  contents. 
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TEMPERATURE  C*C> 

Fig.  1.  Total  particle  concentrations  vs.  temperature 
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Fig.  2a.  Percent  of  particles  <  1  n:ir.  diameter  vs.  temperature 


temperature  cc) 


Fig.  2c.  Percent  of  particles  2-3  nun  diameter  vs. 


temperature 


temperature  Cc> 


Percent  of  particles  nm  diameter  vs.  temperature 


PERCENT  OF  PARTICLES  4 -5  MM  MAM 
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Fig.  2e.  Percent  of  particles  k-5  run  diameter  vs.  temperature 


PERCENT  OF  PARTICLES  >  5  KM  WAM 
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